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Abstract 

We study the dynamics of entanglement for the XY-model, one-dimensional spin systems coupled 
through nearest neighbor exchange interaction and subject to an external time-dependent mag- 
netic field. Using the two-site density matrix, we calculate the time-dependent entanglement of 
formation between nearest neighbor qubits. We investigate the effect of varying the temperature, 
the anisotropy parameter and the external time-dependent magnetic field on the entanglement. 
We have found that the entanglement can be localized between nearest neighbor qubits for certain 
values of the external time-dependent magnetic field. Moreover, as known for the magnetization 
of this model, the entanglement shows nonergodic behavior, it does not approach its equilibrium 
value at the infinite time limit. 
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I. INTRODUCTION 



Quantum entanglement is regarded as the resource of quantum information processing 
with no classical analo g 1 ! 2 ' 3 ' 4 . The corresponding investigation is currently a very active area 
of researc h 5 ' 6 ' 7 ' 8 ' 9 ' 10 due to its potential applications in quantum communication, such as 
quantum teleportatio n 11 ' 12 , superdense coding 13 , quantum key distribution^ 4 -, telecoloningi^ 
and decoherence in quantum computer o 16 ' 17 . 



Multiparticles systems are the central interest in the field of quantum information, in 
particular, the quantification of the entanglement contained in quantum states, because the 
entanglement is the physical resource to perform some of the most important quantum infor- 
mation tasks, like quantum information transfer or quantum computation. Osterloh et.al^ 
connected the theory of critical phenomena with quantum information by exploring the en- 
tangling resources of a system close to quantum critical point in a class of one-dimensional 
magnetic systems. 

Recently^, we have demonstrated that for a class of one-dimensional magnetic systems 
entanglement can be controlled and tuned by varying the anisotropy parameter in the 
XY Hamiltonian and by introducing impurities into the systems in the equilibrium state. 
However, offering a potentially ideal protection against environmentally induced decoher- 
ence is difficult in information encoding and readout. An important motivation to study 
the dynamics of entanglement while varying an external magnetic field is to investigate 
whether it is possible to protect the entanglement from the effects of environment such as 
an external magnetic field and change of temperature. 

Amico et. al2S study the dynamics of entanglement in one-dimensioanl spin systems 
using Ising-type models. They analyze the time evolution of initial Bell states created in 
a fully polarized background and on the ground state. They have found that the pairwise 
entanglement propagates with a velocity proportional to the reduced interaction for all the 
four Bell states. Moreover, they show that the "entanglement wave" evolving from a Bell 
state on the ground state turns out to be very localized in space time. 

In this paper, we consider the dynamics of a set of localized spin- 1/2 particles coupled 
through an exchange interaction and subject to an external time-dependent magnetic field. 
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In Sec. II, we introduce the Liouville equation for the density matrix and present the 
numerical solution of the general XY-model in a lattice with N sites in an external time- 
dependent magnetic field h(t). In Sec. Ill, the solutions presented in the previous section are 
used to compute the magnetization and spin-spin correlation functions. The entanglement 
of formation is briefly introduced in Sec. IV and expressed in terms of the different spin-spin 
correlation functions. Finally, Sec. V is devoted to the results and discussions. 



II. SOLUTION OF THE TIME-DEPENDENT XY MODEL 

In this section, we present the numerical solution of the XY model for a one-dimensional 
lattice with N sites in an external time- dependent magnetic field h(t). The Hamiltonian for 
such a chain of interacting spins, with nearest-neighbor interaction only, is given by 

j N j N N 

H = - j (1 + 7) E « i " 2 " 7) E ' E h ®°i ' W 

i=l ~ i=l i=l 

where J is the coupling constant, a a are the Pauli matrices (a = x, y, z), and 7 is the degree 
of anisotropy. We set J = 1 for convenience. The periodic boundary condition is cyclic, 
namely, a a N+l = of. 

The standard procedure used to solve Eq.Q is to transform the spin operators into 
fermionic operator o 2 ^ 2 7 . Let us define the raising and lowering operators af, Oj, 

4 = \(°? + i<>f), ai=\{of-i<$), (2) 
in terms of which the Pauli matrices are given by 

°-=at + a u of = ^^, of = 2a+a 4 -J. (3) 

These operators can be expressed in terms of Fermi operators b iy bf 

i-l i-l 

a,i = exp(—ni ^ b~j~bj)bi, af = bfexp(ni ^ bfbj). (4) 
j=i 3=1 

Next, we introduce the Fourier transform for a general h(t) 

1 n/2 1 N/2 

bj = -7= exp(ij(j) p )c+ , bj = -= J2 exp(-ij<j) p )c p , (5) 

V iV p= _ N / 2 ViV p= _ N / 2 
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where <p p = Thus, the Hamiltonian assumes the following form 



N/2 

H = ocp{t)[c^c p + c+ p c_ p ] + i5 p {c+ct p + c p c^ p ] + 2h(t) . 
P =i 



(6) 



where a p {t) = —2coscf) p — 2h(t) and 5 P = 2jsincj) p . Since [H p , H q ] = , we can write Eq. 
(0 as 



with 



N/2 

P =i 



H p = a p {t)[c+c p + c+ p c„ p ] + i5 p {c+ct p + c p c_ p ] + 2h(t) 



(7) 



(8) 



This means that the space of H can be decomposed into noninteracting subspaces, each of 
four dimensions. Using the following basis for the pth. subspace: 



(|0>;c+c+ p |0>;c+|0>;c± p |0>), 



(9) 



we can explicitly obtain 



/ 2h{t) 



H P (t) 



-i$v 







° \ 



iS p -4cos<j>p - 2h(t) 
-2cos</> v 

\ -2coscj> p ) 



(10) 



In this paper, the initial condition chosen at t = is thermal equilibrium of the system, 
namely, the density matrix of the pth subspace at time t p p (t) is given by 



p p (0) = e -^(°), p = 1/kT 
k is the Boltzmann constant. Therefor, using Eq. (fTTH) . we obtain 



(11) 



Pp(0) 



2/3cos^>+2/3A[h(0)] 



/ k{ 2 \ 

^21 ^22 ^ " 

fcg 3 

V fcf 4 / 



(12) 



where 



A[/i(0)] = {[cos<f) + h(0)} 2 + fsin 2 ^ 2 , 

and the matrix elements are given by 

{A[/i(0)] + cos0 + /i(0)}e- 4/3A [ /l (°)] + {A[h(0)] - cos<f) - h(0)} 



^ii 



2A[/i(0)] 



(13) 



(14) 
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ft 21 



_^{l_ e -4/3A[MO)]} 



^22 



4A[M0)] 

{A[/i(0)] - cos<p - h{0)}e~ A ^ h( -^ + {A[/i(0)] + cos<P + /i(0)} 



2A[/i(0)] 

fc P = „-20A[fc(O)] 



Let C/ p (t) be the time-evolution matrix in the pth subspace, then 

dU p {t) 



dt 



Since HJt) is in a block diagonal form 



U p (t)H p (t) , ft = 1 



(15) 

(16) 

(17) 
(18) 

(19) 



U p (t) 



(u{ x Uf 2 \ 
U$ 2 
[/ 3 P j 



(20) 



V U^J 

where the upper-left block is determined from 

2h(t) -i8p 
iS p —4cos(p p — 2h(t) 

Thus, the Liouville equation of the system which is given by 

.dp(t) 



( u p u 


u? 2 ) 


( D?i 






u p 22 




u% 2 ) 



(21) 



dt 



[H(t),p(t)} 



(22) 



can be solved exactly because it can be decomposed into uncorrelated subspaces. In the pth 
subspace, the solution of Liouville equation is 

p p (t) = U p (t)p p (0)U p (ty . (23) 

In this study the magnetic field will be presented by a step function of the form, 



h(t) 



a t<0 
b t>0 



which will allow us to obtain the solution of Eq. (fTTJj) 



U p (t) = e 2itcos ^ 



U2i U 22 

u% 3 
V t/f 4 / 



(24) 
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where 



-i(cos</> + b)sin[2tA(b)] 

MP) + 

_ -6sin[2tA{b)} 

U 12 — ' 



cos[2tA(b)} 



TjP 



21 



u. 



22 



2A(6) 

Ssin[2tA(b)} 
2A(b) 

i(coscj) + b)sin[2tA(b)] 

TjV — TJP _ 1 



+ cos [2tA {b)\ , 



(25) 

(26) 

(27) 

(28) 
(29) 



From Eq. (|2~3j) we can get 



p p (t) = e 2/3cos0+2/3A[h(O)] 



where the matrix elements are given by 



P21 P22 
pf 3 

V p^ 4 / 



(30) 



p _ {5 2 (b - a)sin 2 [2tA(b)] + C} e - 4 ^ A ( Q ) + S 2 (a - b)sin 2 [2tA{b)] + rj 
Pl1 ~ 4A 2 (6)A(a) (31) 

5(1 - e- 4 ^){A(b)sin[AtA{b)]{b - a) + i{A 2 {b) + 2(a - &)(cos0 + 6)sm 2 [2tA(6)]}} 



P12 



v 

P22 



4A 2 (6)A(a) 

{5 2 {a - b)sin 2 [2tA{b)] + v } e - 4 ^ + <5 2 (6 - a)sin 2 [2tA(b)} + ( 

4A 2 (6)A(a) 

A. = (p? 2 )* 

Pis = P P u = e" 2/3A(a) 



(32) 
(33) 

(34) 

(35) 



with 



C = 2A 2 (6)(A(a) +cos(p + a) 
7] = 2A 2 (6)(A(a) -cos<f)-a) 



(36) 
(37) 
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III. MAGNETIZATION AND SPIN-SPIN CORRELATION FUNCTIONS 



The magnetization in the XY model is defined as 



which can be written in terms of the operators c p , c p as 



N/2 



M 



N 



P =i 



where M p = c p c p + c± p c_ p — 1 . So we can get the z-direction magnetization 

M (t) = IZM^P] 1 ^ Tr[M pPp (t)} 
ZU N Tr[p] Tr[p p (0)} ' 

Using Eqs. Q,P and 0| Eq. fP gives 

J_ ^ tanfr[/?A(a)]{2<S 2 (6 - a)sin 2 [2tk(b)} + 4A 2 (6)(cos0 + a)} 



P =i 



A 2 (6)A(a) 



The three instantaneous spin-spin correlation functions are defined as 

QX , qX qX . qU , Qj/ rry ^ Q 2 ^ Q 2 Q 2 ^ 

°«m — ^ °l °m ^ i °lm ~ ^ J i D m ^ > — ^ °« °m 

Lieb Schultz and Mattis (LSM)^ show that 



qx 



lm 



qy 
°lm 



- < BiAi + iBi...A m _iB m _iA m > 



-(—1) m < AiBi + iAi + iBi +2 ...B m ^iA m ^iB m > 



Sim — ^ < AiBiA m B m > 



where 



(38) 



(39) 



(40) 

(41) 
(42) 

(43) 
(44) 
(45) 

(46) 



Ai = bf + h- B i = bf - h. 

These three correlation functions are given as expectation values of products of fermion 
operators. Using the Wick 2 ^ theorem, the expressions can be expressed as Pfaffians (pf). 
In particular, we have 



/ < B t A l+1 > < B t B l+1 > 
< A l+1 B l+1 > 



-,Pf 



< BiBm-! > < B L A m > \ 

< Ai +1 B m -i > < A l+1 A m > 



(47) 
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/ < A t B l+1 > < A t A l+1 > 
< B l+1 A l+1 > 



i_r\l-m 
Sf =>— -1 pf 

lm A 



< AiA m -i > < A L B m > \ 
< Bi+iA m —i > < Bi +1 Bm > 



< B m — \A m -\ > < Bra—xBm > 

< A m - X B m > / 



< AiBi > < A t A m > < AiBm > 
;/'./ I /•'•>•• • ■ BiB m > 



(48) 



(49) 



where 



1 ^ 2 1 27r(m — Z)» 

+ 2(a - 6)(cos0 + 6)szn 2 [2tA(6)]} + cos (^^-^) 

{5 2 (6 - a)sm 2 [2iA(6)] + 2A 2 (6)(cos0 + a)}[l - e " 4/3A(a) ]} , (50) 



1 ^(m-ZV i5(a-b)sin(^^)sin[4tA(b)}tanh((3A(aj)^ 

< m >= ivg {2c ° g( iV } + A(&jAW } ' 

(51) 

1 ,27r(m-0*\ ^(a-6)sm(2^^)sm[4tA(6)]tan/i(^A(a)) 1 

< ^ m >= iV^ { - 2c ° s( N } + ApW } • 

(52) 



IV. ENTANGLEMENT OF FORMATION 

The concept of entanglement of formation is related to the amount of entanglement 
needed to prepare the state p, where p is the density matrix. It was shown by Wootters^ 
that 

E(p)=S(C(p)) J (53) 

where the function £ is given by 



£ = K ), (54) 

where h(x) = —xlog2X — (1 — x)logi(\ — x) and the concurrence C is defined as 

C(p) = max{0, Ai — A 2 — A3 — A4 }. (55) 



For a general state of two qubits, Aj's are the eigenvalues, in decreasing order, of the Her- 
mitian matrix 

R= ^VppVp, (56) 

where p is the density matrix and p is the spin-flipped state defined as 



P = {<Jy <g> Oy)p*(a y ® Oy). (57) 

Alternatively, the Aj's are the square roots of the eigenvalues of the non-Hermitian pp. 
Since the density matrix p follows from the symmetry properties of the Hamiltonian, the p 
must be real and symmetrical, plus the global phase flip symmetry of Hamiltonian, which 
implies that [ofcrJ,p] = 0, we obtain 



P 



f Pi,i PlA ^ 

P2,2 P2,3 

P2,3 P3,3 

\Pl,4 p M J 



(58) 



with 



K = v/pi,lP4,4+|pi,4 1, A fo = v /p2,2P3,3+|P2,3|, K = \ VPlAPM ~ |Pl,4 1 1 , A d = 1^2^3,3-^2,31- 

(59) 

Using the definition < A >= Tr(pA), we can express all the matrix elements in the 
density matrix in terms of the different spin-spin correlation functions: 



Plil = ~Mf + X -M* m + S? m + i (60) 

P2, 2 = \m z x - \M z m - Sf m + i, (61) 

Ps,3 = \M z m - \Ml - S't m + \, (62) 

p 4 ,4 = ~Mf - \Ml + Sf m + i (63) 

P2, 3 = Sf m + Sf m} (64) 

Pi,4 = Sf m - Sf m . (65) 
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V. RESULTS AND DISCUSSIONS 



Our goal is to examine the dynamics of entanglement in the presence of varying external 
magnetic field, temperature and the anisotropy parameter 7. First we describe the dynamics 
for the Ising model with 7 = 1. For a constant magnetic field, it is convenient to define 
a dimensionless coupling constant A = J/h. This model is known to undergo a quantum 
phase transition at A c = 1. The magnetization < a x > is different from zero for A > 1 
and it vanishes at the transition^. However, the magnetization < a z > is different from 
zero for any value of A. At the quantum phase transition the correlation length diverges as 
£ ~ |A — A c | _1 . When A — > 0, the ground state becomes a product of spins pointing in the 
positive z-direction. However, in the limit A — > 00, the ground state becomes again a product 
of spins pointing in the positive x-direction. In both limits the ground state approaches a 
product state, thus the entanglement vanishes. When A = 1, a fundamental transition 
in the form of the ground state occurs and the system develops a nonzero magnetization 
< a x >^ which grows as A is increased. The calculations of entanglement show that it 
rises from zero in the two limits A — > and A —>■ 00 to a maximum value near the critical 
point A c = 1. Moreover, the range of entanglement, that is the maximum distance between 
two spins at which the concurrence is different from zero, vanishes unless the two sites are 
at most next-nearest neighbors. 

In Fig (1) we show how the nearest neighbor concurrence C(i,i + 1) evolves with time 
when the initial C(i,i + 1) close to the maximum. We choose the parameters a = b = 1.001, 
a = b = 0.5 and the step function with a = 1.001 and b = 0.5. Thus at t = 0, A close to one 
and C(i,i + 1) close to maximum. As time evolve, C(i,i+ 1) oscillate, but it does not reach 
it is equilibrium value at t — ► 00. Barouch et. al.— have shown the nonergodic behavior 
of the the magnetization for the XY-model. The limit t — > 00 of the magnetization does 
not approach its equilibrium value. This phenomenon, the magnetization is not an ergodic 
observable in this model was discusses earlier by Mazui^. The concurrence C(i, i + shows 
a similar behavior, that is nonergodic, since it is related to the magnetization and spin- 
spin correlation functions. In the lower panel of Fig. (1), we calculate the thermal nearest 
neighbor concurrence C(i,i + 1) as a function of time t for kT = 0.5 and kT = 1.0. For this 
model, the entanglement is nonzero only in a certain region in the (kt — A) planed- The 
entanglement is largest in the vicinity of the critical point A c = 1 and kT = 0, this is the 
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quantum critical regime. As expected the concurrence decreases with increasing temperature 
at A close to one and the oscillations disappeared at kT = 1.0. 

In Fig. (2) we show the nearest neighbor concurrence C(i,i + 1) as a function of time t 
at kT = and kT = 0.5 for different parameters of the magnetic field, a step function with 
an initial field a = 0.5 and a final field b = 5.0. For a = b = 0.5, A = 2 > A c = 1 and for 
a = b = 5.0, A = 0.2 < A c . One can see that the concurrence starts oscillations when the 
external magnetic field is applied and reaches a limit when t — > oo , which is again not the 
equilibrium limit. 

In order to investigate the property of concurrence at equilibrium, we calculate the three 
spin-spin correlation functions as defined in Eq. (J43j) . Eq . (J44j) . Eq. (J45|) and the magnetization 
in Eq. ([38)1 . Figures (3) and (4) show the behavior of three spin-spin correlation functions 
and the magnetization as a function of time t at kT = and kT = 0.5 respectively. As 
reported by Barouch et. al.— , the magnetization of the Ising model does not approach the 
equilibrium state limit. Furthermore, we find that the three spin-spin correlation functions 
do not approach the equilibrium state at t — > oo. 

To show the effect of the initial and final external magnetic field strengths on the 
entanglement with t — > oo, we show in Fig. (5) the nearest neighbor concurrence C(i,i + 1) 
as a function of the parameters a and b at kT = and j — 1. For a < 1 region, the 
concurrence increases very fast near 6=1 and reaches a limit C(i, i + 1) ~ 0.125 when 
b — > oo . It is surprising that the concurrence will not disappear when b increases with 
a < 1. This indicates that the concurrence will not disappear as the final external magnetic 
field increase at infinite time. It shows that this model is not in agreement with the obvious 
physical intuition, since we expect that increasing the external magnetic field will destroy 
the spin-spin correlations functions and make the concurrence vanishes. In our previous 
calculations^, we have found that the concurrence approached a maximum when the 
external magnetic field is close to the critical point. The concurrence approaches maximum 
C(i,i + 1) ~ 0.258 at (a = 1.37,6 = 1.37), and decreases rapidly as a ^ b. This indicates 
that the fluctuation of the external magnetic field near the equilibrium state will rapidly 
destroy the entanglement. However, in the region where a > 2.0, the concurrence is close 
to zero when b < 1.0 and maximum close to 1. Moreover, it disappears in the limit of b — > oo. 
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Recently, it was reported that the nearest neighbor concurrence will decrease as the 
temperature increases^ at the equilibrium state. It is interesting to investigate the effect 
of temperature on the concurrence in our model. Fig. (6) shows the nearest neighbor 
concurrence C(i,i+ 1) as the parameters a and b varies at kT = 1. The concurrence in the 
region where a < 1 disappears, and the sharp peak shown at kT = decrease by increasing 
the temperature. The maximum C(i,i + 1) ~ 0.195 occurs at (a = 1.76, b = 3.0). 

We now move to consider the dynamics for the anisotropy parameter 7 / 1. In Fig. (7) 
we show the nearest neighbor concurrence C(i,i + 1) as a function of time t at kT = 
for the same parameters of the magnetic field shown in Fig. (2) for 7 = 1. One can see 
that the concurrence starts oscillations when the external magnetic field is applied and 
reaches a limit when t — > 00 . As for the case 7 = 1, this limit is not equivalent to the 
concurrence for the equilibrium state with a = b = 5.0. In the lower panel we calcu- 
late the thermal nearest neighbor concurrence C(i, i + 1) as a function of time t for kT = 0.5. 

Up to now we examined the dynamics of nearest neighbor concurrence C(i, i + 1) as a 
function of the magnetic filed parameters (a,b), the temperature and anisotropy parameter 
7. To describe the dynamics of the next nearest neighbor concurrence C(i, i + 2), first we 
compare in Fig. (8) the behavior of C(i,i + 1) and C(i,i + 2) as a function of time for 
same parameters a = b = 1.15 at kT = and 7 = 1. Although C(i,i + 2) shows oscillatory 
behavior as for C(i,i + 1), but the magnitude is very small compared with the C(i, i + 1). 
Moreover, by increasing the temperature C(i,i + 2) decreases and vanishes for kT > 0.125 
as shown in Fig. (9). In Fig. (10) we show the dynamics of C(i,i + 2) for 7 = 0.5 at 
kT = and kT = 0.1. For this case the value of C(i,i + 2) is larger than the case with 
7 = 1 but with similar dynamics. It is interesting to mention that C(i,i + 2) is different 
from zero along the magnetic field parameters a = b as shown in Fig. (11) for 7 = 1. The 
maximum C(i,i + 2) ~ 0.004 occurs at (a = 1.0, b = 1.0). For 7 = 1 we need to consider 
only the dynamics of the nearest neighbor C(i, i + 1) and next nearest neighbor C(i, i + 2) 
since C(i, i + 3) vanishes. 

In summary, we have studied the dynamics of entanglement for one-dimensional spin 
systems in an external magnetic field of a step function form. We observed that the en- 
tanglement shows nonergodic behavior. Due to the coherence of the pairwise entanglement 
with the environment, the change of external magnetic field decreases the nearest and second 
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nearest pairwise entanglement. However, at low temperatures, we have found that there are 
some regions where there is a decoherence of the entanglement due to the change of the 
external magnetic field. Finally, we have found that an increase of the temperature in the 
system will always decrease the pairwise entanglement. 
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FIG. 1: The nearest neighbor concurrence C{i,i + \) for different external magnetic field strengths 
a and b as a function of time t for kT = and kT = 0.5 and 7=1. 
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Fig. 2 
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FIG. 2: The nearest neighbor concurrence C{i,i + \) for different external magnetic field strengths 
a and b as a function of time t for kT = and kT = 0.5 and 7 = 1. 
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Fig. 3 




FIG. 3: The spin-spin correlation functions and the average magnetization per spin for different 
external magnetic field strengths a and b as a function of time t at kT = for 7 = 1. 
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Fig. 4 
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FIG. 4: The spin-spin correlation functions and the average magnetization per spin for different 
external magnetic field strengths a and b as a function of time t at kT = 0.5 for 7 = 1. 
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FIG. 5: The nearest neighbor concurrence C(i,i + 1) as functions of different external magnetic 
field strengths a and b at time t — ► oo and temperature kT = for 7 = 1. 
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Fig. 6 




FIG. 6: The nearest neighbor concurrence C(i,i + 1) as functions of different external magnetic 
field strengths a and b at time t — ► oo and temperature kT = 1.0 for 7 = 1. 
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Fig. 7 
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FIG. 7: The nearest neighbor concurrence C{i,i + \) for different external magnetic field strengths 
a and b as a function of time t for kT = and kT = 0.5 and 7 = 0.5. 
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Fig. 8 
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FIG. 8: Comparison of the nearest neighbor concurrence C(i,i + 1) and the next nearest neighbor 
concurrence C(i,i + 2) for different external magnetic field strengths a and b as a function of time 
t for kT = and 7 = 1. 



23 



Fig. 9 
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FIG. 9: The next nearest neighbor concurrence C(i,i + 2) for different external magnetic field 
strengths a and 6 as a function of time t for kT = 0.1 and kT = 0.125 for the case 7=1. 
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Fig. 10 
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FIG. 10: The next nearest neighbor concurrence C(i,i + 2) for different external magnetic field 
strengths a and b as a function of time t for kT = and kT = 0.1 for the case 7 = 0.5. 
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Fig. 11 




FIG. 11: The next nearest neighbor concurrence C(i,i + 2) as functions of different external 
magnetic field strengths a and b at time t — ► oo and temperature A;T = for 7 = 1. 
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